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We discuss the statistical properties of lifetimes of electromagnetic eigenmodes in dielectric microresonators 
with fully chaotic ray dynamics. Using the example of a resonator of stadium geometry, we find that a recently 
proposed random-matrix model very well describes the lifetime statistics of long-lived resonances, provided 
that two effective parameters are appropriately renormalized. This renormalization is linked to the formation 
of anomalously short-lived resonances, a mechanism also known from the fractal Weyl law and the resonance 
trapping phenomenon. 
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I. INTRODUCTION 

Over the past decade optical microresonators have attracted 
considerable attention in various fields of physics [1]. Mi- 
croresonators are basic components for single-photon emit- 
ters H, ultralow threshold lasers JSllj], and solid-state cav- 
ity quantum electrodynamics la 0], to name just a few 
of the possible applications. Deformed dielectric microdisk 
resonators are of particular interest since they allow for di- 
rectional light emission, an effect that can be understood and 
optimized by explorin g fu ndamental principles of ray-wave 
correspondence fi H, [HI [H O Ej- This directionality 
can be employed to design microlasers which make best use 
of the total output intensity. Whether a carefully designed 
mode admits substantial lasing depends on how its lifetime 
compares to the lifetimes of other modes in the amplifica- 
tion window. This information is encoded in the imaginary 
parts of the resonance frequencies. However, the associated 
spectral statistics have been studied only for a few special 
geometries, such as the circular disk without and with sur- 
face roughness ITil [T3l. while other investigations focussed 
on disordered medi a Ilia, \Y7n an d resonators with small open- 
ings H0 m EH 

In order to address the spectral properties of dielectric mi- 
crodisk resonators from a broader perspective it is desirable to 
adopt concepts which have been successfully applied to other 
classes of open systems. In particular, the characterization 
of spectral statistics is often based on two ingredients: (a) A 
Weyl law M25I1 . which concerns the average number of modes 
within a given frequency interval, and (b) random-matrix the- 
ory (RMT) [26], which addresses the statistical fluctuations 
and correlations of these modes. Traditionally both concepts 
have been applied to systems with small openings, for which 
the Weyl law is identical to that of the closed counterpart, 
while RMT can be justified because modes become strongly 
mixed on the scale of the typical life time. 

A framework how to adapt these concepts to dielectric mi- 



crodisk resonators, who are relatively more open systems al- 
lowing for refractive escape along the entire boundary, has 
only emerged very recently. Due to the short typical life- 
time, strong mode mixing is only observed for a reduced 
number of modes. For systems with ballistic escape, where 
the mode-mixing time scale is given by the Ehrenfest time 
||27[ l28i l29l l30ll and the support of the long-living states is 
the fractal repeller, it has been established that the number of 
modes that exceed a certain lifetime is given by a fractal Weyl 
law HI mHHH HI. Based on a generalization of the 
repeller for refractive escape, the significance of the fractal 
Weyl law to dielectric microdisk systems with small refrac- 
tive index has been demonstrated for the stadium geometry in 
Ref. ll36ll . In the present work we consider the same geometry 
and explore the predictive power of RMT, whose extension to 
dielectric microresonators was proposed in Ref. ll37ll . The ap- 
plication of this extension requires to determine two effective 
parameters (a time scale for the internal dynamics and an ef- 
fective number of modes), which are sensitive to the formation 
of short living resonances (a mechanism also at heart of the 
fractal Weyl law lO . as well as the resonance trapping phe- 
nomenon Bal ). We concentrate on moderately large values 
of the refractive index, for which the short-living resonance 
states resemble the bouncing-ball modes of the closed stadium 
billiard [39], and find that the numerical lifetime statistics of 
long-living modes agree very well with the RMT predictions. 

The outline of this paper is as follows. The background 
Section|II]provides a brief introduction into the spectral prop- 
erties of dielectric microdisk resonators and the dynamical 
model of Ref. 1I37I1 . In Section [III] we present numerical re- 
sults for complex resonance frequencies in stadium-shaped 
microresonators and then analyze the statistics of long-living 
resonances in two steps. First we show that the mean proper- 
ties of the resonance frequencies can be characterized in terms 
of the two effective parameters mentioned above. Based on 
these parameters, we secondly show that fluctuations of the 
life times are well captured by RMT. The results are summa- 
rized in the concluding Sec. [IV] 
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II. THEORETICAL BACKGROUND 
A. Dielectric microdisk resonators 

Dielectric microdisk resonators effectively have a 2D ge- 
ometry, defined by a region of constant refractive index n 
which is embedded into vacuum (where the refractive index is 
unity). The classical dynamics inside such a resonator is sim- 
ilar to a billiard, where a point-like particle moves freely in a 
two-dimensional domain with elastic reflections at the bound- 
ary 114011 . Depending on the shape of the boundary a billiard 
can show a variety of dynamical behaviors ranging from inte- 
grable to fully chaotic 114 ill . Light rays in a microdisk behave 
similarly as they are totally reflected at the interfaces as long 
as the angle of incidence is larger than the critical angle for 
total internal reflection. If, however, the angle of incidence 
is smaller than the critical angle, light can escape refractively 
according to Snell's and Fresnel's laws. Hence, dielectric mi- 
crodisk resonators are leaking billiards JUIHH^]. 

The leakiness of the resonator limits the photon lifetime, 
which can be described by associating complex values to the 
resonance frequencies u>. These frequencies can be obtained 
from the Helmholtz equation [43] 



(1) 



where c is the speed of light in vacuum. For a piecewise con- 
stant refractive index, this equation holds for both transverse 
magnetic (TM) and transverse electric (TE) polarization. For 
TM polarization the electric field is perpendicular to the res- 
onatorplane with E z = Re[ip(x, y)e~ iuJt ], and the wave func- 
tion tp and its normal derivative d±tp are continuous across the 
boundary of the resonator. For TE polarization, ip represents 
the z-component H z of the magnetic field, and the functions 
ip and n~ 2 d±vp are continuous across the boundary. For plane 
interfaces, these boundary conditions yield the Fresnel reflec- 
tion amplitudes 
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for resonator photons with dimensionless transverse momen- 
tum p = sin x, where \ is tne angle of incidence. 

At infinity, outgoing wave conditions are imposed which 
results in quasi-bound states with frequencies to situated in the 
lower half of the complex plane. Whereas the real part is the 
usual frequency, the imaginary part is related to the lifetime 
r = — l/[2 Imu;] . The quality factor of a quasi-bound state is 
defined by Q = — Re ujj [2 Im u] . 

An alternative method to obtain the quasibound states sets 
out with the scattering matrix S, which relates the amplitudes 
of the incoming wave components to the amplitudes of the 
outgoing wave components (see, e.g., Refs. lfl8l [l9l l2ll I22I 
l23ll ). At real frequencies, the scattering matrix is unitary. The 
condition of a quasibound state corresponds to the situation 



where outgoing wave components are admitted in absence of 
any incoming radiation. Under these conditions the scatter- 
ing matrix diverges. The resonance frequencies therefore cor- 
respond to the poles of the scattering matrix in the complex 
frequency plane. 



B. Dynamical model of microresonators 

Random-matrix theory delivers a statistical description of 
complex quantum systems on the basis of just a few system 
parameters. In particular, a classification according to fun- 
damental symmetries such as time-reversal symmetry or spin- 
rotation symmetry allows to identify ten classes, among which 
three are the traditional Wigner-Dyson classes 112611 . while a 
further seven classes apply to systems with chiral or particle- 
hole symmetries 14411 . For each symmetry class one can con- 
sider two basis variants of RMT, which either concerns her- 
mitian matrices H (representing the Hamiltonian) or unitary 
matrices U (representing the scattering matrix or the time evo- 
lution operator). In the absence of magneto-optical effects 
photonic systems are represented by the orthogonal symmetry 
class, in which the hermitian and unitary matrices are sym- 
metric under transposition (H = H T , U = U T ). 

Traditional random-matrix theory does not assume any fur- 
ther knowledge of the system, which leads to a uniquely de- 
fined invariant probability measure (the Haar measure). Phys- 
ical systems, however, often contain additional commonali- 
ties. Dielectric microresonators are a good example, since the 
escape is governed by well structured Fresnel laws of reflec- 
tion and refraction. These laws not only affect the far-field 
radiation characteristics, but also determine the confinement, 
and therefore impact on the spectral properties. It is therefore 
desirable to incorporate the specific confinement characteris- 
tics of dielectric microresonators into RMT. In order to study 
spectral statistics, one furthermore needs to accommodate fre- 
quency correlations associated to the internal dynamics in the 
resonator. Both requirements are met by a variant of RMT 
recently proposed in Ref. [37]. In this model the scattering 
matrix takes the form 
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are determined from the Fresnel reflection amplitudes iff]) at 
discrete impact parameter pi = (21 — M — 1)/M (with 
M=dimS and I = 1, . . . , M) JH. The frequency depen- 
dence is characterized by a time scale r which represents 
the typical time of flight through the resonator. The M- 
dimensional matrix U can be interpreted as a unitary time evo- 
lution operator which propagates photons between encounters 
with the interface. 

In the dynamical model d3j, the poles of the scattering ma- 
trix are determined by the eigenvalue equation 



KUip = exp(-i6)ip, 



(4) 
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FIG. 1: Sketch of a stadium-shaped dielectric microresonator with 
L = R. 



where 6 = lot is a dimensionless quantity. The operator on 
the left is subunitary, and therefore each resonance frequency 
in general possesses a negative imaginary part related to the 
lifetime of the state. Eigenvalue problems of this type with 
various specifications of TZ have recently attracted consider- 
able attention (H [B |H 53 IH IU . 

In the RMT version of the dynamical model Q, U is taken 
from the circular orthogonal ensembles. The statistical de- 
scription is then completed by the specification of the two pa- 
rameters M and r. Once these parameters are determined, one 
can investigate statistical fluctuations. We now implement this 
strategy for a specific model system. 



III. LIFE TIME STATISTICS IN THE STADIUM 
MICRORESONATOR 

A microdisk resonator with stadium-shaped cross-section 
is illustrated in Fig. Q] The region with refractive index n is 
bounded by two semicircles with radius R and two straight 
lines of length L. The corresponding closed cavity, the sta- 
dium billiard, is a paradigm for fully developed chaos j50l,l5lll . 
We consider the case L = R, for which the closed classical 
billiard has the largest Lyapunov exponent [52]. 

The open stadium resonator has been studied both theo- 
retically and experimentally in a number of works, with fo- 
cus on the spatial mode structure and ray-wave correspon- 
dence ll53l l54l l55i l56l l57ll . In a first step towards the charac- 
terization of the resonance frequency spectrum of the stadium 
resonator, the applicability of the Weyl law has been investi- 
gated in Ref. 13611 . 

We start our statistical analysis of the spectral fluctuations 
with the numerical determination the resonance frequencies. 
In the next step, we relate the averages of the real and imag- 
inary parts of these frequencies to the two parameters r and 
M of the dynamical model ([3]). This fixes the RMT version 
of this model and delivers a prediction for the distribution of 
life times, which we compare in the final step to the actual life 
time statistics of the stadium resonator. 
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FIG. 2: Position of dimensionless complex resonance frequencies 
with TM polarization for the stadium microresonator with refractive 
index n = 5. 



A. Resonance frequencies 

The resonance states in the stadium billiard can be classi- 
fied by their parity with respect to the horizontal and vertical 
symmetry lines of the stadium billiard. In the computations 
this can be exploited by considering four desymmetrized ver- 
sions of the resonator, which differ by the boundary conditions 
imposed on the two symmetry lines (Dirichlet or Neumann 
boundary conditions for states of odd or even parity with re- 
spect to the given symmetry line, respectively). As explained 
in detail in Ref. [ 36] the boundary element method 15811 and 
the harmonic inversion technique 115911 are very efficient tools 
to compute the quasibound modes. The combination of these 
two methods allows us to gather a large number of resonant 
frequencies to, as is required for a reliable statistical analysis. 
A dimensionless form of these frequencies can be obtained 
by considering the scaled frequencies £1 — luR/c. For the 
present study we collected resonance frequencies in the inter- 
val 66 < n Re Q < 82.5 for both polarizations (TM and TE), 
all parities, and two different values of the refractive index, 
n = 3.3 (for GaAs) and n = 5. The latter value we have cho- 
sen since a broad interval of refractive indices allows a more 
convincing validation of our theory. 

For illustration, the position of the dimensionless resonance 
frequencies in the complex plane for n = 5 and TM polar- 
ization (all parities) is shown in Fig. [2] The figure shows 
a clear separation of short-living resonances which are orga- 
nized along well-structured bands (Im SI < —0.015), and long 
living resonances whose position does not suggest any clear 
organization principle (ImO > —0.015), and therefore call 
for a statistical analysis. 

FigureOa) and (b) shows examples of long-lived and short- 
lived modes in the chaotic microstadium. The complex spa- 
tial pattern of the long-lived mode is in strong contrast to 
the rather simple pattern of the short-lived mode. Short-lived 
modes in the microstadium resemble bouncing-ball states in 
the closed stadium billiard B39I1 . Bouncing -ball modes are 
usually ignored in the statistical analysis of the quantum spec- 
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FIG. 3: (Color online). Near-field pattern of a long-lived (a) and a 
short-lived resonance mode (b) with Q = 15.71886 — iO. 00337 and 
Q, = 15.74567 — z0.0378. As in Fig. [2] these modes are computed 
for TM polarization and n — 5. 



trum of closed systems ll6()ll . In the same spirit we restrict 
in the following our attention to the long-living resonances, 
which are selected via a cut-off of the decay rate (ImO > 
-0.04 for n = 3.3, TM; ImO > -0.1 for n = 3.3, TE; 
ImO > -0.015 for n = 5, TM; ImO > -0.025 for n = 5, 
TE). 



B. Determination of characteristic parameters 

A useful reference point for the quantification of the spec- 
tral fluctuations of the long-living resonances are the fre- 
quency and time scales set by their mean spacing along the 
real frequency axis, and their mean decay rate. The nu- 
merical values of the dimensionless quantities (Af2)i tot - > and 
— (ImO) s are collected in Tab. [I] Here the subscript V refers 
to the stadium billiard, while the superscript 'tot' for the spac- 
ing indicates that we do not discriminate the resonances by 
parity. We now describe how this information can be used 
to determine the parameters M and r which serve as input 
of the RMT model. Because of parity conservation, we are 
interested in the effective parameters per symmetry class. 

Let us first assume that we could ignore the existence of 
well organized sequences of short-living resonances, which 
corresponds to the case of a fully wave-chaotic resonator. The 
mean resonance spacing then follows the ordinary Weyl law, 



(AO) 



Wcyl 



2tt R 2 
i 2 Ren A' 



(5) 



where A = (4 + ir)R 2 is the area of the billiard. Analogously, 
the dimension M( tot ' of the scattering matrix can be estimated 



TABLE I: The first two rows compare the dimensionless mean fre- 
quency spacing of long-living resonances in the stadium resonator in 
the interval 66 < Re(nf2) < 82.5 (irrespective of parity) to the pre- 
diction from the Weyl law {5}. The second and third row tabulate the 
mean decay rate of the long-living resonances and the mean dimen- 
sionless decay rate in the dynamical model (fJJ (obtained in RMT in 
the asymptotic limit M S> 1). As described in the text, the relative 
scaling of these quantities can be used to determine the time scale 
r [Eq. l(9j] and effective matrix dimension M (for fixed parity) [Eq. 
lIlOH, which are tabulated in the last two rows. 
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by applying the Weyl law to the perimeter S 
the billiard, 



AL 



(tot) 
Weyl 



nReO S 
n R' 



(27T + 4)i?0f 



(6) 



In the dynamical model, eigenvalue equation has M so- 
lutions with uniform statistical distribution of the phase Re 9. 
This translates into a mean dimensionless frequency spacing 



(AO) 



(tot) _ 



2tt R 



(7) 



M( tot ) TC 

The combination of Eqs. (0, ©, and results in the relation 



nix A 



T , 



(8) 



which coincides with Sabine's law for the mean propagation 
time between encounters with the boundary. If applied to 
resonances of a fixed parity, the frequency spacing AO = 
4(AO)( tot ) quadruples, while the effective number of modes 
M = M ( tot ) /4 is reduced by a factor four. It follows from Eq. 
(0 that the time scale r is not affected by these substitutions. 

The finiteness of the time to establish wave chaos in an 
open system results in the reduction of the number of long- 
living quasibound states, which is compensated by the forma- 
tion of short-lived states 13311 . A striking example is the emer- 
gence of a fractal Wey l law for systems with a finite Ehren- 
festtime JM HI Hi HI HI- The formation of short-lived 
states also often results in resonance trapping, a phenomenon 
which reduces the average decay rate of the long-living reso- 
nances 1 380 . The existence of well-organized bands of reso- 
nances deep in the complex-frequency plane (see Fig. |2]l in- 
dicates that both mechanisms are at work in the stadium res- 
onator. As we now show, the consequences can be described 
by renormalized parameters M and r. 

Since the reduced number of long-living resonances and 
resonance trapping both impact on the average characteristics 
of complex frequencies (affecting the real and imaginary part, 
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n=3.3, TM 



stadium 
RMT 




0.08 



0.08 



FIG. 4: (Color online). Distribution of decay rates (in dimensionless 
units) for resonances in the stadium resonator with refractive index 
n = 3.3, compared to the prediction of RMT based on the dynamical 
model {3} with T — 1.44ro and M — 38. Short-living resonances 
in the stadium which form bands in the complex-f2 plane (see Fig. 
|2j are neglected via a cut-off. Top panel: TM polarization, Im Q > 
—0.04. Bottom panel: TE polarization, Im > —0.1. The large 
ticks indicate the value of the average dimensionless decay rate (see 
Tab. IB. 



respectively), the renormalized parameters can be determined 
without resorting to the statistical fluctuations. As a first rela- 
tion, we compare the mean decay rate of long-living states in 
the stadium to the prediction of RMT, 



CT 



(lm9) 



RMT 



(ImQ) s 



(9) 



For the second relation we evaluate Eq. (0 using the nu- 
merically observed frequency spacing of the long-living reso- 
nances, M s (tot) = [27r/(An)i tot) }[R/cT]. The effective num- 
ber of modes per symmetry class of fixed parity is therefore 
given by 



M 



R 



2(Afi) 



(tot) 



(10) 



In principle, Eqs. (0 and ( TTOb need to be solved self- 
consistently because the RMT average (Im S)rmt also de- 
pends on M. For large values of M, however, the RMT av- 
erage saturates. Practically, therefore, t can be obtained from 
Eq. ® by fixing M to a large value. The actual value of M 
then follows from (TTOb . This procedure is consistent if the re- 
sulting value M 1. Table|T]tabulates the values for r and M 



n=5, TM 
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RMT 




0.02 
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FIG. 5: (Color online). Same as Fig. [4] but for higher refractive 
index, n = 5. The cut-offs are ImS) > -0.015 (TM) and Imfi > 
-0.025 (TE). 



which result when this procedure is applied to the resonances 
of the stadium resonator. 

For all four combinations of refractive index and polariza- 
tion, we find that the collision time r differs from the predic- 
tion To of a totally wave-chaotic resonator by a factor of order 
unity, while the effective number of modes M is also almost 
identical in all cases. Taking the mean over all four cases we 
obtain the values 



t = 1.44 t , M — 38. 



C. Life time statistics 



(ID 



When a resonator is employed as a microlaser, the modes 
with the smallest decay rate within the amplification window 
tend to be the first that cross the lasing threshold. In a simple 
model, lasing occurs as soon as the gain counters the losses 
through the interfaces iflrjl [l9l l2ll l22l l23ll . The analysis of the 
mean decay rate in the preceding section therefore delivers a 
rough estimate of the necessary gain at the lasing threshold. 
The distribution of resonances in the complex plane (Fig. [2]), 
however, suggests that the decay rates fluctuate broadly. Such 
large fluctuations entail that some modes have a much smaller 
decay rate, which results in a reduced lasing threshold. We 
now turn to the question whether these fluctuations are cap- 
tured by the RMT variant of the dynamical model ||3). 

Histograms of the decay rates of the resonances in the sta- 
dium resonators are shown in Figs. [4] and [3] For comparison, 
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the figures also show the prediction of the RMT model, with 
parameters given by Eq. ( fTTT i. Since these parameters apply 
to the quasibound states of a fixed parity, the total spectral 
statistics are obtained by the superposition of four systems 
with identical effective parameters, but independent micro- 
scopic realizations of the Af-dimensional random matrix U. 
The RMT results in Figs.|4]and|5]are based on 50000 realiza- 
tions per combination of polarization and refractive index. 

In all cases we find excellent agreement of all key features, 
such as the modal value of the decay rate, the width of the 
peak around the modal value, and the shape of the distribu- 
tion in the two tails. In particular, it should be noted that the 
average decay rates (tabulated in Tab. UJ and shown as large 
ticks in the figures) are noticeably larger than the modal value 
of the decay-rate distribution. This difference arises from the 
long tail of resonances with relatively large decay rate, which 
closely follows a slow power law. Such power-law tails are a 
robust feature of RMT models [ 18], including models formed 
on the basis of the eigenvalue problem (0) [48]. Further- 
more, in the RMT model the width of the distribution around 
the modal value depends noticeably on M. The good agree- 
ment of the width with the actual histograms lends direct sup- 
port to the procedure employed in the determination of the 
renormalized effective parameters ( fTTT i. We emphasis that the 
renormalized parameters are determined independently and 
do not result from a fitting to the numerically computed life- 
time statistics. 



IV. SUMMARY 

Dynamical models with scattering matrices of the form 
(0 have recently attracted considerable attention, both within 
RMT [46l|47l,|48l] as well as for individual dynamical systems. 
In particular, these models have been explored in association 
with the Ehrenfest time 11281 l29l 13011 and the related question 
of fractal Weyl laws Il33[l34[r35ll49ll . Some investigations of 



the Ehrenfest time employ an effective RMT model for the 
scattering matrix irJll l62l l63ll . In these investigations it was 
found that certain transport properties (such as universal con- 
ductance fluctuations and shot noise) are faithfully described 
by effective RMT, while others (such as the weak-localization 
correction) are not l29ll . 

The present work is motivated by the suggestion in 
Ref. B37I1 that the dynamical model (0 can also be used to 
describe the spectral fluctuations of long-living resonances in 
chaotic dielectric microresonators. We combined this model 
with a variant of effective RMT, which is directly applied to 
the internal round-trip operator U l33ll . 

We confirm that the RMT variant of the dynamical model 
captures the key features of the life time statistics in a spe- 
cific dielectric microresonator, the stadium-shaped resonator. 
This comparison requires a careful determination of the ef- 
fective collision time r of photons with the interfaces and the 
effective number of modes M that are mixed by the reflec- 
tions at the boundary. Our analysis shows that these param- 
eters can be reliably obtained by using the mean resonance 
spacing and decay rate of the long-living resonances. Com- 
pared to the assumption of complete wave chaos, the effective 
parameters are renormalized in a uniform way. This renor- 
malization finds its natural explanation in the formation of 
anomalously short-lived resonances that resemble bouncing- 
ball modes, which changes the mean resonance spacing of 
long-living resonances IB3I1 . and the resonance trapping phe- 
nomenon [38], which changes the mean decay rate. 
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